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Abstract 



We investigate universal features of the off-equilibrium sequential and conservative 
fragmentation processes with the dissipative effects which are simulated by the Gaus- 
sian random inactivation process. The relation between the fragment multiplicity 
scaling law and the fragment size distribution is studied and a dependence of scal- 
ing exponents on the parameters of fragmentation and inactivation rate functions is 
established. 
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Fragmentation is the universal process which can be found at all scales in the nature. 
The most general sequential binary and conservative fragmentation processes with scale- 
invariant fragmentation and inactivation rate functions, have been previously studied in 
much details • The phase diagramme of these off-equilibrium processes has been 

established and the universal aspects of both the fragment size distribution and the total 
number of fragments distribution (i.e., the multiplicity distribution) have been determined 
H|4[] . In this fragmentation-inactivation binary (FIB) model Jl],|2| , one deals with fragments 
characterized by some conserved scalar quantity that is called the fragment mass . The 
anscestor fragment of mass N is fragmenting via an ordered and irreversible sequence of 
steps. The first step is either a binary fragmentation, (N) — > (J) + (N — j) , or an 
inactivation (N) — > (N)* . Once inactive, the cluster cannot be reactivated anymore. The 
fragmentation leads to two fragments, with the mass partition probability ~ Fj,N~j ■ I n the 
following steps, the process continues independently for each active descendant fragment 
until either the low mass cutoff for further indivisible particles (monomers) is reached or all 
fragments are inactive. For any event, the fragmentation and inactivation occur with the 
probabilities per unit of time ~ Fj^-j and ~ Ik respectively. The fragmenting system and 
its evolution is completely specified by that rate functions and the initial state. It is also 
useful to consider the fragmentation probability pf without specifying masses of descendants 
: PF(k) = J2i=i Fi,k-i {Ik + YnZi Fi^-i)' 1 ■ If the instability of smaller fragments is smaller 
than instability of larger fragments, PF{k) is an increasing function of fragment mass and 
the total mass is converted into finite size fragments. This is the shattered phase. The 
fragment mass independence of Pf(^) at any stage of the process until the cutoff-scale for 
monomers characterizes the critical transition region. The multiplicity anomalous dimension 
: 7 = d(ln < m >)/d(lnN) , is the order parameter in the FIB model. It equals 1 in the 
shattering phase and takes the intermediate value between and 1 in the critical transition 
region. 

For most fragmenting systems, the off-equilibrium relaxation process ceases due to a 
dissipation. The dissipation is not always scale-invariant as considered in Ref. 1 but, on 
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the contrary, it is often characterized by a definite and usually small length scale. It is 
then an open question to which extent the fragmentation processes which on one side are 
driven by the homogeneous scale-invariant fragmentation rate function and on other side 
are inactivated at a certain fixed scale by the random inactivation process, may develop 
scale-invariant and universal features in both the fragment mass distribution n(k) and 
the fragment multiplicity distribution P(m) . This question is important in view of the 
widespread occurence of scale-invariant fragment mass distributions n(k) ~ k~ T and the 
lack of convincing arguments for using homogeneous dissipation functions in many processes 
including parton cascading in the perturbative quantum chromodynamics (PQCD) || or 
the fragmentation of highly excited atomic nuclei, atomic clusters or polymers. In this 
work, we address this fundamental question using the FIB process with the homogeneous 
fragmentation rate function : Fj,k-j = [j{k — j)] a , and with the dissipation at small scales 
which is modelled by the Gaussian inactivation rate function : 

h = cexp[~ ] . (1) 

An asymptotic (t — > oo) fragment mass distribution in the critical transition region 
of FIB model with scale-invariant dissipation phenomena , is a power law with an 
exponent r < 2 . In the shattering phase, the fragment mass distribution is also power law 
but with an exponent r > 2 . Another characteristic observable is the fragment multiplicity 
distribution : P{m) = J2k -Pfc( m ) , where Pfe(m) is the probability distribution of the number 
of fragments of mass k . This quantity has been intensely studied in the strong interaction 
physics . Of particular importance is a possibility of asymptotic scaling of multiplicity 
probability distributions : 

< m > s P(m) = , z {5) = - - — s (2) 

where the asymptotic behaviour is defined as<m>-^oo,m^oo for a fixed 

(m/ < m >) - ratio. < m > is the multiplicity of fragments averaged over an ensemble 

of events. The scaling law @ means that for example data for differing energies (hence 
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differing < m >) should fall on the same curve when < m > 5 P(m) is plotted against the 
scaled variable z/$\ = (m— < m >)/ < m > s . Some time ago Koba, Nielsen and Olesen 
(KNO) suggested an asymptotic scaling (fj) with 5 = 1 in the strong interaction physics 
|7]] . The same scaling has been found also in the critical transition region of scale-invariant 
FIB process for pp > 1/2 and a > -1 [|. Recently, Botet, Ploszajczak and Latora (BPL) 
reported another scaling limit in (Q) with 5 = 1/2, which holds in the percolation and in 
the shattering phase of scale-invariant FIB process 0] . 5 = 1/2 and 1 are the two limiting 
values since 8 > 1 or <5 < 1/2 are incompatible with the scaling hypothesis (§) . 

The study presented in this Letter correspond to the domain a > — 1 of fragmentation 
rate functions Fj,k-j . Many known homogeneous fragmentation kernels correpond to this 
domain. These include the singular kernel a = — 1 in the PQCD gluodynamics || , 
a = —2/3 for the spinodal volume instabilities in three dimensions Q , a = +1 in the 
scalar A0g field theory in six dimensions , and many others . For a < — 1 , the 
fragmentation process is dominated by the splitting (k) — » (k — 1) + (1) at each step in the 
cascade, and leads to the finite limiting value of < m > independently of the initial size 
iV . In this evaporation phase, the scaling solution (g) does not hold and the multiplicity 
anomalous dimension is equal zero when iV — > oo . This phase is not relevant for the problem 
we want to address in this Letter. 

Without restricting the generality of our discussion, we will present below results for 
fragmentation kernels with : a = — 1 and a = +1 . The upper part of Fig. 1 shows 
multiplicity distributions for a = — 1 in the scaling variables (0) for 5 = 1 (the upper 
left part), and fragment mass distributions for the same parameters (the upper right part). 
The cascade equations of Gaussian FIB model have been solved by Monte-Carlo simulations 
|1|J^] for different initial system sizes (N = 1024, 4096) and for the following exemplaric 
parameters : c = 1 and a = 0.1 , 1 of inactivation rate function = /^(c, cr) . We have 
made exhaustive analysis of P(m) for a broad range of c, a parameters , finding in all cases 
the KNO scaling (5 = 1). We have found the KNO scaling uniquely for a = — 1 . The shape 
of KNO scaling function ^(z^) depends on the precise value of both c and o . 
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In the lower left part of Fig. 1, we show typical multiplicity distributions for a = 
+ 1 which are plotted for different system sizes in the BPL scaling variables (5 = 1/2 ). 
The corresponding fragment mass distributions are shown in the lower right part of Fig. 
1. Again, the precise form of BPL scaling function $(z(i/2)) depends on the chosen set of 
parameters c and a . In contrast to these results of Gaussian FIB model, fragmentation 
process in the scale-invariant FIB model for any value of exponent a may be found either 
in the critical transition region or in the shattering phase depending on the homogeneity 
index (3 of the inactivation rate function Ik = I\hP [|l]|| . This means that e.g. both for 
a = —1 and +1, one may see either the KNO scaling or the BPL scaling of multiplicity 
distributions depending on the precise value of the homogeneity index of the inactivation 
term. 

Concerning the fragment mass distributions, Fig. 1 shows the distributions for a = 
— 1,+1 and different parameters of Gaussian inactivation rate function Ik{c, a) . For o 
larger than ~ 0.5 , one finds the power law distribution of fragment masses for any value of 
parameter c . In the studied case : a = 1, c = 1 , the exponent r equals 1.8 and 2.8 for 
a = — 1 and a = +1 respectively. For a given a , the value of exponent r is remarkably 
independent of a but depends strongly on the value of parameter c in I k (c,a) . For a 
smaller value of o {a — 0.1 is shown in Fig. 1), the fragment mass distribution decreases 
exponentially and the shape of scaling function resembles the Gaussian distribution. The 
form of this exponential distribution depends both on c and a parameters. 

As a generic case for a = — 1 , we have found the scale- invariant region of power law 
fragment mass distributions with r < 2 for o above ~ 0.5 , and the exponential region of 
mass distributions for a less than ~ 0.5 . The power law region is completely analogous to 
the critical transition region of scale-invariant FIB model for a > — 1 and pp > 1/2 Jl|-§[ , 
because the multiplicity anomalous dimension in both models is : 



We have verified validity of this relation in Gaussian FIB model for a broad range of c, a val- 



7 = t — 1 




(3) 
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ues. In the exponential region , 7 is always equal 1 independently of the value of parameter 
c , i.e. this region is in the shattering phase. One should remind that shattering in the scale- 
invariant FIB model is related exclusively with the BPL scaling , whereas in the Gaussian 
FIB model for a = — 1 the KNO scaling holds. 

The fragment size distributions for a = +1 and different values of a behave similarly 
as for the a = — 1 case, except that now for a above ~ 0.5 the power law exponent is 
t > 2 . For all a, i.e. in both exponential and power law regions of mass distribution, 
the multiplicity anomalous dimension is 7 = 1 and the BPL scaling holds. This generic 
situation is completely analogous to the multiplicity behaviour found in the shattering phase 
of scale-invariant FIB model . 

Whenever the fragment size distribution is a power law, the KNO scaling of multiplicity 
distributions is associated with r < 2 and the BPL scaling of multiplicity distributions 
with r > 2 in both scale-invariant and scale-dependent regimes of dissipation. This clearly 
indicates a direct relation between the multiplicity scaling law and the fragment mass distri- 
bution scaling regimes in the FIB model. In view of the generality of FIB process, it would 
be very interesting to test this relation experimentally. A novel aspect of the Gaussian FIB 
model is associated with properties of multiplicity scaling in the new region of exponential 
fragment mass distributions. In this region, BPL scaling holds for a = +1 whereas KNO 
scaling is seen for a = —1 . 

In Fig. 2 we plot for different values of the parameter c the normalized cumulant factorial 
moment of order two |10[ : 72 = (< m{m — 1) > — < m > 2 )/ < m > 2 , vs the width a of 



inactivation rate function Ik{c, a) . The exponent of homogeneous fragmentation kernel is : 
a = —1 . For this choice of a , the KNO scaling holds and 72 becomes the second moment 
of scaling function which is independent of initial mass iV |6]|| . For each point (c , 

cr) , cascade equations of the FIB model have been solved exactly by recurrent formula || 
up to initial system size N = 2 18 . As can be seen in Fig. 2, the multiplicity fluctuations as 
measured by 72 are extremely small in the exponential region for a less than ~ 0.5 . The 
change of 72 when passing from the power law to exponential region is continuous but the 
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largest variations of 72 (er) appear at o ~ 0.5 . For large values of a , the cumulant factorial 
moment approaches a limiting value which depends on the value of parameter c . 

The experimental informations about 72 are not numerous and concern mainly charged 
particle multiplicities at relativistic and ultrarelativistic energies. The DELPHI Collabora- 
tion reported the data on hadron production in e + e~ annihilations for the center of mass 



(cm.) energy of y/s = 91GeV finding 72 = 0.04 JTTJ . In hadron-hadron collisions 7r + — p 



K + —p,p — p,p — p for cm. energies ranging up to 1000 GeV f||l2ll , values of 72 increase 



from about 0.05 to 0.3 as energies increase to collider values. Distribution of galaxy counts 
in the regions of sky covered by the Zwicky catalogue |l3j yields 72 — 0.3 |14|]. Indepen- 
dently of the question whether the KNO scaling holds in all those different physical systems 
, the measured values of 72 clearly exclude the exponential region of Gaussian FIB process. 
Much more information could be extracted if in addition to the moments of the multiplicity 
distribution also the mass distribution would be available. In high energy lepton and/or 
hadron collisions for example , this would require measuring the hadron mass distribution. 

In conclusion, we have demonstrated that the off-equilibrium binary fragmentation with 
scale-invariant fragmentation kernel and the scale-dependent inactivation simulating the 
dissipation at small scales, yields the fragment mass and fragment multiplicity distributions 
which are scale-invariant for a broad range of parameters. This is an important finding 
because most of fragmentation processes in nature which have these scale-invariant features 
are probably not associated with the dissipative processes acting at all scales. The scale- 
dependent fragmentation processes may also develop strong scale-invariant fluctuations (the 
KNO scaling) though the region of their appearance is restricted to the particular value 
of exponent : a = — 1 , of the homogeneous fragmentation function. The region at a = 
— 1 and a above ~ 0.5 is the critical transition region of Gaussian FIB process. For 
other values of a , the fragment multiplicity distributions obey the BPL scaling, i.e. the 
small amplitude limit of scaling multiplicity fluctuations. Another transition zone of the 
Gaussian FIB model is defined by the width a of inactivation rate function. At a ~ 0.5 , 
the fragment size distribution changes from exponential (for o < 0.5) into power law (for 



a > 0.5) . The form of scaling function $>(z$) , together with the form of fragment mass 
distribution n(k) impose strong constraints on the choice of basic functions of FIB kinetic 
equations : the fragmentation and inactivation functions. This has been demonstrated on 
the example of hadron production data in the e + e~ annihilation [H| . Results of this Letter 
show that the closing of gap between experimental observables related to the fragment mass 
distribution and/or the fragment multiplicity distribution and the basic ingredients of the 
kinetic theory, i.e. the rates of activation Fj^-j and inactivation , can be achieved for 
many physical systems in the nature. 
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Figure captions 
Fig. 1 

Multiplicity probability distributions in the scaling variables (see eq. (0)), and the fragment 
mass distribution for two homogeneous fragmentation kernels and two Gaussian inactivation 
rate functions. Each set of data corresponds to 10 6 independent events of Monte-Carlo 
simulations. 

(i) Upper left part : the fragmentation kernel with a = — 1 and the inactivation rate 
function (|l|) for c = 1 and two typical values of a. Two sets of data are plotted for two 
different total mass : N = 1024 (crosses) and N = 4096 (circles). These data are plotted in 
the KNO form, i.e. : 5 = 1 (see eq. (0)). 

(ii) Upper right part : the fragment mass distributions in a double-logarithmic scale are 
shown for the same parameters a,c,a as in (i). The total mass is N = 4096. Big stars 
represent results obtained for the same value of a, c parameters and for a much larger 
value of cr (a — 10 ) , to show the independence of the scaling part of the fragment mass 
distribution with the value of a . The line in between points is shown to guide the eyes. 

(iii) Lower left part : the same as in (i) but for the fragmentation kernel with a = +1 . 
These data are plotted in the BPL form, i.e. : 5 = 1/2 (see eq. @). 

(iv) Lower right part : the fragment mass distributions for a = +1. Parameters c,a,N as 
in (ii). 

Fig. 2 

The cumulant factorial moment 72 of the fragment multiplicity distribution is plotted vs 
the width parameter a of the Gaussian inactivation function ([!]) with c = 0.5,1,5 . The 
homogeneous fragmentation kernel is taken with a = —1 . Each point corresponds to system 
of size N = 2 18 , and the values of 72 are calculated by solving exact recurrent equations. 
The line joining points is shown to guide the eyes. 
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